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Sofiane M. Boucenna, 1,2 ’* Jose W.F. Valle, 1, 1 and Avelino Vicente 1,3 ’ t 

1 Instituto de Fisica Corpuscular (CSIC-Universitat de Valencia), Apdo. 22085, E-46071 Valencia, Spain. 
2 INFN, Laboratori Nazionali di Frascati, C.P. 13, 1000)4 Frascati, Italy. 

:l IFPA, Dep. AGO, Universite Liege, Bat B5, Sart-Tilman B-4000, Liege 1, Belgium 

(Dated: September 10, 2015) 

The simplest realization of the inverse seesaw mechanism in a SU(3)c 0 SU(3)l ® I/(l)x gauge 
theory offers striking flavor correlations between rare charged lepton flavor violating decays and the 
measured neutrino oscillations parameters. The predictions follow from the gauge structure itself 
without the need for any flavor symmetry. Such tight complementarity between charged lepton 
flavor violation and neutrino oscillations renders the scenario strictly testable. 
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Preliminaries 

Beyond the discovery of the Higgs boson [1, 2] no 
signs of genuine new physics have shown up so far at 
high energies. However, the existence of new physics has 
been established with the discovery of neutrino oscilla¬ 
tions [3, 4], implying the existence of lepton flavor viola¬ 
tion and nonzero neutrino masses. Unraveling the origin 
of the latter constitutes one of the main challenges of par¬ 
ticle physics. While the prevailing view is that neutrino 
masses arise from physics associated with unification, 
they might as well signal novel TeV-scale physics lead¬ 
ing to potentially large charged lepton flavor violating 
(LFV) rates and possibly also new phenomena testable 
at the LHC [5]. In this case it could well be that new 
physics may actually show up mainly in the form of lep¬ 
ton flavor violation, boosting the motivation to search for 
charged LFV phenomena such as the rare decay /i — > eq. 
In fact, the current limit BR(/z —i eq) < 5.7 x 10 -13 [6] 
already puts severe constraints on models of new physics. 

Models based on the SU (3)c 0 SU (3) l 0 U( 1).y gauge 
theory (3-3-1) constitute a minimal extension of the 
Standard Model (SM) that accounts for the existence 
of three families of fermions, the same as the number 
of colors [7, 8]. They provide an economical scheme 
to generate tiny neutrino masses radiatively from TeV 
scale physics [9] and could lead to successful gauge cou¬ 
pling unification through neutrino masses and TeV scale 
physics [10]. Moreover, they naturally solve the strong 
CV problem by including in an elegant way the Peccei- 
Quinn symmetry [11, 12]. 

Here we focus on the phenomenology of lepton flavor 
violation in the SU (3 )c 0 SU (3)i0(7(l)x schemes. For 
definiteness we focus on the simplest implementation of 
the inverse seesaw mechanism within the 3-3-1 model. 
We show that it offers striking flavor correlations between 
rare charged lepton flavor violating decays and the mea¬ 
sured neutrino oscillations parameters. Such predictions 
result from the gauge theory structure itself without the 
need for imposing any specific flavor symmetry. We ana¬ 
lyze the complementarity between charged LFV and neu¬ 


trino oscillations, a feature which may render the 3-3-1 
scenario strictly testable within the upcoming generation 
of LFV searches. 


The Model 

We consider a variant of the model introduced in [9] 
in which neutrinos get masses via the inverse seesaw 
mechanism instead of quantum corrections. The model 
is based on the SU (3)c 0 SU (3)l 01/(1)x gauge symme¬ 
try, extended with a global U{\)c so as to consistently 
define lepton number 1 . We also invoke an auxiliary 
parity symmetry in order to ensure a realistic quark 
mass spectrum. The model contains three generations 
of lepton triplets (V’l), two generations of quark triplets 
(Ql’ 2 )> one generation of quark anti-triplet (Q£), along, 
of course, with their iso-singlet right-handed partners, 
and accompanied by three generations of neutral fermion 
singlets (S). The gauge symmetry breaking is imple¬ 
mented through three scalar anti-triplets (<(> 1 , 2 , 3 ) ■ The 
particle content of the model is summarized in table (I). 
The fundamental fermions interact through the exchange 
of 17 gauge bosons: the 8 gluons of SU(3)c, the 8 “weak” 
W, bosons associated to SU(3)l (4 of which form 2 
electrically charged bosons, and the rest are neutral), 
and the B boson associated to U(l)x- 

The lepton representations in table (I) can be decom¬ 
posed as: 

/ £- \ e ^’ T 

i>L = ( I , (1) 

\ NC )l 

where we identify 7V£ = (vr) c [8]. 


1 For other inverse seesaw constructions within 3-3-1 scenarios see 

[13, 14], 
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TABLE I: Particle content of the model. Here Ur = 
{ur, cr, tfl}, Dr = {dR, SR,bR} and = ( d' R ,SR ). 


In the scalar sector, on the other hand, we have: 



( 2 ) 


After electroweak symmetry breaking, the electric charge 
and lepton number assignments of the particles of the 
model follow from the action of the operators: 

Q = T 3 + Dr + X ; (3) 

i = 75 n + c - (4) 

The relevant terms in the Lagrangian for leptons are: 
-£lep = y e $LlR4>l+y a 4 , L^L4>l+y S $LS(l)2+ 1 -^ S C S+ h.c. 

(5) 

where y l and y s are generic 3x3 matrices, while y a is 
anti-symmetric and ms is the 3x3 Majorana mass term 
for the singlets S (symmetric, due to the Pauli principle). 

Scalar potential and symmetry breaking 

The scalar potential of the model can be written as: 

v = +a,; | 4 +y> J -i0ii%i 2 

* 

+f (010203 + h.c.) + ml (0 2 03 + h.c.), (6) 

where /ti j2 , 3 , f and m s are parameters with dimensions 
of mass. The two latter couplings break the Z 2 softly. 
For simplicity we denote all the dimensionless couplings 
by A and take m s = 0. 

In full generality, the scalars of the model are al¬ 
lowed to take vacuum expectation values (VEVs) in 
the following directions (0i) T = (fci, 0,0)/\/2, (0 2 ) T = 
(0,fc 3 ,n)A/2, an d (0 3 ) T = (0 ,k,2,n')/\/2. However, in 
order to recover the SM as a low energy limit, we as¬ 
sume the hierarchy ^ 1 , 2,3 n,n'. Moreover, we assume: 

&3 = n' = 0, which together with the Z 2 symmetry guar¬ 
antees the existence of a simple pattern of realistic quark 
masses (see below). 


We define the covariant derivative in the usual way as: 

D l j = d lx -i J2 > ( 7 ) 

groups 

where is the gauge boson, T a are the generators of 
the group and the sum extends over all gauge groups 
included in SU( 3) c (g) SU(3)l <g> U(l)x- Assuming k\ ~ 
L ’2 = k <C n, and keeping only the leading order terms, 
one finds that the mass spectrum of the charged scalars 
is given as: 2 


M 2 (0±) 

= 0 , 

(8) 

M 2 ((0± + 0±)/ v / 2) 

= 0, 

(9) 

M 2 ((0±-0±)/ v / 2) 

1 , 

~ V2 fn ' 

(10) 

M 2 (0±) 

~ y/2fn. 

(11) 


On the other hand the masses of the neutral CP-even 
scalars are, up to corrections of 0(k 2 ): 


Af 2 (5i(0? + 0°)/v / 2) 

~ (2X + V2-- 

n 

I i -)k 2 

(12) 

Af 2 (3?(0? - <fP 3 )/y/ 2) 

1 , 

~ t n , 

V2 

(13) 

M 2 (SR0°) 

= 0 , 

(14) 

M 2 (SR0 2 ) 

~ 2A n 2 , 

(15) 

M 2 (3?0 3 ) 

~ V2 f n. 

(16) 


Finally, the masses of the neutral CP-odd scalars at lead¬ 
ing order are given as: 


Af 2 (9(0° + 0°)/v / 2) 

~ V2fn, 

(17) 

Af 2 (3(0°-0°)/y/2) 

= 0 , 

(18) 

M 2 (S0°) 

= 0 , 

(19) 

M 2 (90 2 ) 

= 0, 

(20) 

M 2 (90 3 ) 

~ ~^=fn. 

V2 

(21) 


The massless scalars found in the above equations corre¬ 
spond to the degrees of freedom ‘eaten-up’ by the charged 
and neutral gauge bosons, respectively, which acquire the 


2 We identify the corresponding (approximate) eigenstates be¬ 
tween parentheses. 
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following masses: 


m^y 

1 2 7 2 

= 2 9 ^ ’ 

(22) 

TOpp/ 

1 2 2 
= - a 92U , 

(23) 

to| 

ff 2 2 (4«?? + 3 5 2 ) 2 

2 (<7i + 3g 2 ) ’ 

(24) 

m Z' 

= \{9l + ?>9l)n 2 , 

(25) 

m'x 

= m Y =\ 9l n 2 ■ 

(26) 


Notice that since <f is singlet under the SU{2) R sub¬ 
group contained in SU(3)l, the VEV n will control 
the four new gauge bosons masses and break SU (3) l to 
SU(2)l■ On the other hand, SU(2) R ® U(l)y is broken 
at the electroweak scale by the k\ and fc 2 VEVs down 
to the electromagnetic U(1)q symmetry. For / ~ n all 
the scalars of the model are naturally heavy, except one 
state that we can identify with the SM Higgs boson, i.e., 
H = (</>? + 4$)/V2, in good approximation. Indeed, its 
couplings to the fermions confirm that the state H is the 
one that gives mass to SM fermions. 

Quark sector 


Neutrino masses and inverse seesaw mechanism 

The presence of the small term S C S, in eq. (5), ex¬ 
plicitly breaks U(l)c and provides the seed for lepton 
number violation leading to neutrino masses via the in¬ 
verse seesaw mechanism. Indeed, after spontaneous sym¬ 
metry breaking of the electroweak gauge group, we get 
the following 9x9 neutrino mass matrix, in the basis 
(v,N,S) [15]: 

/ 0 m D 0 \ 

M = 0 Ml, (30) 

V m s) 

where = '/2k\y a and M = ~^ny s 3 - The in¬ 
verse seesaw-induced light neutrino masses can be writ¬ 
ten as [15]: 

m„ = mo (M T ) 1 ms M -1 ■ (31) 

Here, the matrix M can be taken diagonal without loss 
of generality. Using this freedom and taking into account 
that my) is anti-symmetric, eq. (31) can be expressed in 
terms of an effective symmetric 3x3 matrix, M _1 
\! ms M -1 , as: 

m v = —mo M _1 ms M _1 mo = — mo M _1 mo ■ (32) 


We now turn to the quark sector. From the symme¬ 
tries of the model, see table (I) , it follows that the quark 
Lagrangian is given by: 


A simple implication of the antisymmetry of the “Dirac” 
entry my, is that Det(m.„) = 0, so that the lightest neu¬ 
trino in this model must be massless at the tree level. 


Quarks = Q 1 i 2 y u U R cj ) * l +Q 1 L 2 y d D R( t > l + Q 1 L 2 y d d R ^ 2 
+ Q\ y u U R (f>3 + Q\ y d D R (f >i + Q! y u t' R 4>2 
+ h.c., (27) 


where we defined d R = ( d' Rl s ' R ), U R = {u R , c Rl t R } and 
D r = {d R: s R ,b R }. This Lagrangian leads to the follow¬ 
ing mass matrices: 


M d = - 


1 

71 


( 2/11 ^2 

2/12 ^2 

2/13 

0 

0 \ 

2/21 

922 k 2 

2/23 k 2 

0 

0 

2/ii *i 

$12 k l 

$13 k l 

0 

0 

0 

0 

0 

$14 n 

$15 n 

V o 

0 

0 

$24 n 

$25 n J 



( 2/ii ^1 

2/12 fc l 

2/13 fc i 

0 ^ 

1 

2/21 fc l 

922 k l 

2/23 

0 

71 

$11 k 2 

$12 k 2 

$13 k 2 

0 


\ o 

0 

0 

$14™ / 


(29) 


Lepton flavor violation predictions 

Let us now proceed to a simple parameter counting. 
On the left-hand side of eq. (31) one has 5 independent 
complex parameters, since Det(?n„) = 0. In contrast, 
on the right-hand side of eq. (31) one has 9 indepen¬ 
dent complex parameters: 3 in mo, and 6 elements in 
M. Therefore, we have 4 (complex) relations among the 
parameters ( y s ,y a , and ms)- One can choose as free pa¬ 
rameters the 3 off-diagonal entries of M -1 , together with 
a global scaling factor to defined through: 

TOj^ = 77U 1 (to^TO 2 * — TO^TOy- 3 ) . (33) 

From eq. (31), we can see that to scales as y/m^ms/M 2 , 
so that for ms ~ 10 eV, M « 1 TeV, and neutrino masses 
of 0(0.1) eV, we obtain m ss 10 -22 GeV. In contrast, the 


Thanks to the Z 2 symmetry, the SM and exotic sub¬ 
sectors are independent of each other and can be adjusted 
individually to easily obtain a realistic quark sector and 
heavy exotic quarks at the same time. 


3 Note that the matrix in eq. (30) does not depend on the condi¬ 
tions imposed on the VEVs. Indeed, even if k% yl 0 the resulting 
linear seesaw term [16—18] would give only a subleading contri¬ 
bution ~ m v (Myv /n) 2 
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diagonal entries of M~ x are functions of its off-diagonal 
elements and mp- We emphasize that eq. (33) is not an 
ansatz , but the most general solution of eq. (32). 

Such a parameterization makes explicit the direct re¬ 
lation between charged LFV observables and neutrino 
oscillation parameters, which is a characteristic feature 
of our model. Indeed, LFV in this model arises from the 
term: 


- £lfv = y a 4 > IC 1 ip L 4 ‘ i + h.c., (34) 


which depends solely upon the coupling y a , hence m£>. 
Using eq. (33) together with m v = U* m() lag Ul , where 
U v is the leptonic mixing matrix in its standard parame¬ 
terization in terms of three mixing angles and the Dirac 
phase (6), one obtains the relevant coupling for LFV as: 


y 


a 


( 0 yf 2 y ? 3 \ 

[ ~Vl2 0 ^23 I 

\ ~Vl3 ~V23 0 / 


\J A "4m 

\/2k\ fh 


t/Ki (NH) 

V A " 4 m + A ™sol ( IH ) 


(35) 


for normal (NH) and inverse hierarchies (IH). Here the 
parameters are functions of the lepton mixing matrix 
parameters, i.e. = y“ ( 6 l i 2 , 0 i 3 , $ 23 , <5), and are given in 
the Appendix. It is remarkable that the Yukawas y a rel¬ 
evant for determination of LFV rates are, up to a global 
scaling factor, fully determined by the parameters mea¬ 
sured in neutrino oscillation experiments. This allows us 
to make definite predictions for LFV observables that can 
be used to provide an unambiguous test of the model, as 
we show below. 


Radiative A —> £ :i 7 decays 

In order to show the predictive power of the model 
here we focus, for definiteness, on flavour-changing lep¬ 
tonic (radiative) decays. This probe constitutes one of 
the most important tests of new physics and has been 
actively sought after in many experiments. The branch¬ 
ing ratio (BR) of the decay of the charged lepton £) —£ £j 7 
is given as: 


BR(£, £ i7 ) 


ml I (: y a Fy a ) ij I 2 


(36) 


where Ti k is the total decay width of £*,, and F is a 
function that depends on the masses and mixings of 
all the particles running inside the loop (summation 
over the different contributions is implicit here). We 
have three different classes of contributions: i) loops 
mediated by the new heavy gauge bosons. These are 
suppressed due to the large scale of the breaking of 



0.2 0.4 0.6 0.8 

sill 2 013 

FIG. 1: The branching ratio of the decay y —> e 7 ver¬ 

sus sin 0i3 for M = 500 GeV and 1 TeV. We take fh = 
2 x ICG 23 GeV, / = 2 TeV, and n = 3 TeV. The vertical 
band is the 3<r range reported in [22], The other mixing an¬ 
gles are taken within their 3<r range [22]. We also show in the 
lower-right corner the y —> ey branching ratio as a function 
of M, using best-fit values for the mixing angles, as given in 
[ 22 ]. 


SU(3)l compared to ii) contributions from the 

exchange of a charged scalar whose mass is ~ \J f n/y/2\ 
and finally Hi) the “standard” loop, mediated by the SM 
W boson and neutrinos. The latter two contributions 
dominate the £* — > £j 7 amplitude, with relative sizes 
depending on the ratio f/M. 

This branching ratio depends on the neutrino mixing 
parameters, the global scaling factor fh, and on the 
neutrino mass hierarchy. As can be seen in eq. (35), the 
off-diagonal entries in y a are larger in the case of IH 
by a factor ~ sj Arn 2 tm /Am 2 ol with respect to NH so 
that, for the same input parameters, one has larger LFV 
effects in IH. 

We compute the various relevant LFV observables us¬ 
ing FlavorKit [19] 4 . The branching ratio of the decay 
y — » e 7 is shown in figure ( 1 ) as a function of sin $13 
for two different values of the (quasi-Dirac [ 8 ]) right- 
handed neutrino mass M = 500 GeV and 1 TeV. The 
vertical band is the 3cr range reported in [22], whereas 
the other mixing angles are randomly taken within their 
3cr range [22]. This figure has been obtained by varying 
M (by taking different values for the y s Yukawa cou¬ 
plings) for the fixed parameters n = 3 TeV, / = 2 TeV, 
and m = 2 x 10 -23 GeV. We also consider degenerate 
right-handed neutrinos, normal hierarchy for the light 


4 This is a computer tool based on SARAH [20] and SPheno [21], that 
increases their capability to handle flavor observables. 
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FIG. 2: Ratio of the BR of y and r decays to ey, i.e., 
BR.(/r —> e 7 )/BR(r —» ey) versus sin 612 . The vertical band 
is the 3a range given in [22]. The tilted band is obtained by 
varying the other mixing angles within their 3a range [22]. 


neutrinos and a vanishing Dirac CP violating phase. 
One notices that the branching ratio is lower for the 
M = 500 GeV case. This is caused by a partial cancel¬ 
lation between the standard loop, mediated by the W 
boson, and the contribution induced by the exchange 
of charged scalars. This cancellation takes place for 
M ~ 400 GeV and is explicitly illustrated in the lower- 
right corner of figure ( 1 ), where the /i —> ey branching 
ratio is shown as a function of M, using best-fit values for 
the mixing angles [22]. The main message from figure (1) 
is that y —>• ey may take place with sizeable rates, close 
to the current limit, or even larger. Given the expected 
sensitivities of upcoming experiments one finds that the 
detection of this and other muon number violating pro¬ 
cesses might become feasible. 

Since the BR depends on a global multiplicative factor, 
it is interesting to consider the ratio of branching ratios 
of LFV lepton decays. It follows from eq. (36) that: 

BR (4 -7 ijy ) _ ml | ( y a Fy a ) lJ \ 2 / V u 

BR(4->47) rnl\(y“Fy°) kn \yT ek - [ ’ 

For the simplest case of nearly degenerate right-handed 
neutrinos, the F functions are all equal and cancel out 
in the fraction. In this case, eq. (37) depends exclusively 
on the ratios of \y a y a \, i.e., only on the neutrino mixing 
angles. The main advantage of considering the ratio of 
branching ratios and a quasi degenerate spectrum is that 
this leads to clean predictions which do not depend on 
the neutrino mass hierarchy nor the loop functions. In¬ 
deed, in this simplified scenario, by combining eq. (35) 
with eq. (37) we obtain the following predictions: 

BR(p, -7 ey) 

BR(r —> ey) 

BR(r —> ey) 

BR(t —>■ yy) 


Ij/23^13 1 2 

K/12P23 | 2 o 

\y a M 2 ~ ’ 


10 , 


(38) 

(39) 


where we have used the best-fit values for the neutrino 
parameters as derived in the global fit of neutrino 
oscillations given in [22] and set (5 = 0. So, when 
the right-handed neutrino spectrum is degenerate, the 
model predicts BR(y —>• ey) BR(r -7 £,7). Therefore, 
given the expected sensitivities for r LFV decays 5 * * , 
the simple observation of r —► ^7 in one (or several) 
of the near future experiments would rule out our 
simplest degenerate right-handed neutrino hypothesis. 
The viable alternative scenario in such cases would be a 
hierarchical right-handed neutrinos spectrum, impliying 
a non-vanishing contribution of the F loop functions 
in the ratio of BRs. In this case the F functions for 
different flavor transitions can take very different values, 
and thus the ratios in Eq. (39) can clearly depart from 
their predictions in the degenerate scenario. 

As an illustration, in figure (2) we show the ra¬ 
tio of the BR of r and y decays to ey, namely 
BR(^t —7 ey)/BR{r -7 ey) as a function of the solar 
mixing parameter sin $ 12 - The other oscillation param¬ 
eters are varied randomly within their 3 ct ranges [22]. 
Similarly, the ratio of the BR of leptonic r decays, i.e., 
BR(r —7 yy)/BR(r —7 ey) depends mainly on the solar 
mixing parameter. 

For other LFV processes such as ii -7 3 £j and y — e 
conversion in nuclei, our results are qualitatively simi¬ 
lar to the ones found in standard low-scale seesaw mod¬ 
els [24]. Loops including neutrinos give the most impor¬ 
tant contributions, leading to LFV rates comparable to 
the ones for the radiative decay £i -7 tjy. This will be of 
special relevance due to the expected sensitivities in the 
coming experiments [25]. The complete study of all LFV 
processes is, however, beyond the scope of this paper. 


Conclusions and discussion 

In summary, we have shown how a simple extended 
SU(3)c <8> SU(3)l <8> U(1)x electroweak gauge symme¬ 
try implementing the inverse seesaw mechanism implies 
striking flavor correlations between rare charged lep¬ 
ton flavor violating decays and the measured neutrino 
oscillations parameters. The predictions follow simply 
from the enlarged gauge structure without any imposed 
flavor symmetry. Such tight complementarity between 
charged LFV and neutrino oscillations renders the sce¬ 
nario strictly testable. A more detailed study of other 
LFV processes will be taken up elsewhere. The scheme 


5 The expected Belle II sensitivities for r radiative decays are 

around 10 —9 [23], whereas the current MEG bound on BR(// -4 

ey) is many orders of magnitude stronger, BR < 5.7 X ICG 13 . 
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also has a non-trivial structure in the quarks sector 
since, thanks to the anomaly cancelation requirements, 
the Glashow-Iliopoulos-Maiani mechanism breaks down, 
leading to a plethora of flavor-changing neutral currents 
in the quark sector [26, 27]. Last but not least, the model 
presents a rich structure of new physics at the TeV scale 
that could be potentially studied in the coming run of 
the LHC. 
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Appendix: y a Yukawa couplings 

We present in this appendix the expressions for the 
y a Yukawa couplings. Using the definitions of the y a 
elements in eq. (35), we find 


Vi 2 = - (e li cos 6 >i 2 sin 6 »i 3 cos 023 - sin 6 »i 2 sin 6 * 23 ) 2 , (40) 

Vi 3 = e* 5 sin 012 cos 0 i 2 sin 0i 3 cos ( 20 2 3 ) — sin 0 2 3 cos 023 (sin 2 0 i2 — e 2lS cos 2 0i2 sin 2 0i 3 ) , (41) 

i /23 = cos 0i2 cos 0i3 (e I<5 cos 0i2 sin 0i3 cos 023 — sin 0i2 sin 023 ) . (42) 
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